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Abstract 

The electromagnetic analog of an elastic spring-mass network is constructed. 
These electromagnetic circuits offer the promise of manipulating electromagnetic 
fields in new ways, and linear electrical circuits correspond to a subclass of them. 
They consist of thin triangular magnetic components joined at the edges by cylin- 
drical dielectric components. (There are also dual electromagnetic circuits con- 
sisting of thin triangular dielectric components joined at the edges by cylindrical 
magnetic components.) Some of the edges can be terminal edges to which electric 
fields are applied. The response is measured in terms of the real or virtual free 
currents that are associated with the terminal edges. The relation between the 
terminal electric fields and the terminal free currents is governed by a symmetric 
complex matrix W. In the case where all the terminal edges are disjoint, and the 
frequency is fixed, a complete characterization is obtained of all possible response 
matrices W both in the lossless and lossy cases. This is done by introducing a 
subclass of electromagnetic circuits, called electromagnetic ladder networks. It is 
shown that an electromagnetic ladder network, structured as a cubic network, can 
have a macroscopic electromagnetic continuum response which is non-Maxwellian, 
and novel. 
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1 Introduction 

In this paper we introduce a new type of electrical circuit, called an electromagnetic 
circuit, which has the potential at a fixed frequency for providing new and easily ana- 
lyzable ways of manipulating electromagnetic fields beyond those provided by electrical 
circuits, photonic circuits, optical lenses, waveguides, photonic crystals, and transfor- 
mation optics. To construct the electromagnetic circuit we draw upon analogs between 



electromagnetism, elastodynamics and acoustics. Analogs between electromagnetism and 



elastodynamics have a long history [see, e.g. Silva (2007) and references therein] and can 



be understood from the viewpoint of differential geometry (lOziewicz 1994p . It is easy 
to see the connection from the underlying partial differential equations when they are 
written in a form which emphasizes the similarity. Considering for simplicity a locally 
isotropic medium, Maxwell's equations at fixed frequency uj take the form 

D = eE, B = /iH, VxE = icjB, VxH=j-iu;D, (1.1) 

where D(x), E(x), B(x), H(x), and j(x) are the complex electric displacement field, 
electric field, magnetic induction, magnetic field, and free current (the physical fields 
are the real parts of e~*'^*D, e~*'^*E, e~*'^*B, e~*'^*H, and e~*'^*j where t is time) and 
e(x, tu) and /x(x, a;) are the complex electric permittivity and complex magnetic per- 
meability. Here the free current j(x) may represent a single frequency component of a 
time varying ion beam current, or a time varying current generated by an electrochem- 
ical potential. It does not include conduction currents aE, where a is the conductivity, 
that instead are included in the term D = eE, through the imaginary part of £:(x, cj). 
(It is difficult, if not impossible, to distinguish oscillating displacement currents from 
oscillating conduction currents.) The elimination of D, B and H leads to the form 



(Milton, Briane, and Willis 2006) 



d f dE 

Cpqrs^-^ ] + iujjg = -Ul^eEq, (1.2) 



dxp \ dx 
where 

and Cpqm = 1 (—1) if {p,<l,^) is an even (odd) permutation of (1,2,3) and is zero other- 
wise]. This is clearly similar to the form of the equations of continuum elastodynamics 

Cpqrs^r^] + fq = PUq, (1.4) 



dxp \ dx. 

in which u(x) and f (x) are the complex displacement field, and body force (the physical 
fields are u = (e~*'^*u)' and f = (e~"^*f)' where the prime denotes the real part) and 
C(x, cij) is the complex elasticity tensor (incorporating viscosity terms through its imag- 
inary part) and p(x, uo) is the density [which, when it is the effective density tensor of an 
isotropic composite material can be complex and has the same properties as a function 



of a; as e(x, a;): see Milton and Willis (2007) and references therein.] At low frequencies 



one often has the approximation that C(x, a;) ^ C'(x) — iuju{ii) where C'(x) is the real 
component of the elasticity tensor, and ^'(x) is the viscosity tensor, incorporating both 
bulk and shear viscosities. Then, upon introducing the velocity v = du/dt, (11.41) reduces 
to 

+ (1.5) 



dXq'^'^'-'dXr '"'dx/ " 9t 
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Figure 1: Sketch of a two-terminal discrete elastic network, where the terminal nodes 
are represented by the black circles, and the three internal nodes are represented by the 
white circles. In the idealized model these nodes have some mass and are infinitely stiff. 
They are are attached by connecting rods (which act as springs) having no mass and non- 
zero stiffness. Ideally the nodes should have infinitely small diameter and the connecting 
rods should be infinitely thin, although in practice one then has to worry about buckling 
which is a non-linear effect. 



which may be more familiar to readers acquainted with the Kelvin- Voigt model of vis- 
coelasticity. 

As these analogies have been known for a long time it is rather amazing that no elec- 
tromagnetic analog of a spring network with masses at the nodes has been constructed. 
Our electromagnetic circuits are this analog. In elastic networks, as modeled by the 
continuum construction of figure [H the density is concentrated at the nodes, while the 
elasticity is concentrated along the edges. Everything is surrounded by void with C = 
and p = 0. 

Besides the interest of electromagnetic circuits for providing a new way of manipu- 
lating electromagnetic fields there is also a fundamental reason for studying them. It 
is becoming increasingly clear that the usual continuum equations of physics do not 
apply to composite materials built from high contrast constituents and having exotic 
microstructures. Thus, one can obtain materials with macroscopic non-Ohmic, possibly 
non-local, conducting behavior, even though they conform to Ohm's law at the microscale 
(IKhruslov 1978UBriane 1998HBriane and Mazhak 19981 iBriane 20021 |Camar-Eddine and Seppecher 2002 



Cherednichenko, Smyshlyaev, and Zhikov 2006), materials with a macroscopic higher or- 



der gradient or non-local elastic response even though they are governed by usual linear 
elasticity equations at the microscale ([Pideri and Seppecher 1997} IBouchitte and Bellieud 20"02l 
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Alibert, deirisola, and Seppecher 2003^ |Camar-Eddine and Seppecher 2003|) , materials with 



macroscopic behavior outside that of continuum elastodynamics even though they are 
governed by continuum elastodynamics at the microscale (IMilton 2007^ . and materials 



with non-Maxwellian macroscopic electromagnetic behavior ( Shin, Shen, and Fan 2007 ), 

even though they conform to Maxwell's equations at the microscale [see also Dubovik, Martsenyuk, and Sc 

where other non-Maxwellian macroscopic equations are proposed]. 

One would really like to be able to characterize the possible macroscopic continuum 
equations that govern the behavior of materials, including materials with exotic mi- 
crostructures. A program for doing this was developed by Camar-Eddine and Seppecher 
(I2002^ 120031) . and successfully applied to the conductivity and elastic equations in the 
three-dimensional static case, assuming the macroscopic behavior was governed by a sin- 
gle potential or displacement field. The program, in essence, consists of four steps: first 
to show that discrete networks can be modeled by a continuum construction; second to 
characterize all possible responses of discrete networks allowing for part of the network 
to be hidden; third to find the possible continuum limits of these discrete networks; and 
fourth to show that these possible continuum behaviors is all that there can be, even 
when one allows for other, non-network based, microstructures. For the dynamic case, 
at fixed frequency, a complete characterization has been obtained of the possible re- 
sponse matrices of multiterminal electrical, acoustic, and elastodynamic networks, both 
in two and three dimensions, thus meeting the second goal of the program in these cases 



(Milton and Seppecher 2008). For Maxwell's equations at fixed frequency our electro- 
magnetic circuits accomplish, in a formal way, the first goal of the program, and we also 
make progress towards the second goal. At this time it is unclear if our electromagnetic 
circuits are sufficiently rich in construction that their continuum limits can model the 
macroscopic behavior of all other, non-network based, microstructures, and in particular 
the question remains open as to whether the Maxwell equations themselves can be recov- 
ered as a continuum behavior of our electromagnetic circuits. It seems clear, however, 
that many non-Maxwellian continuum behaviors can be achieved (see the concluding 
paragraph of the paper). 

We emphasize that, besides similarities, there are also important differences between 
Maxwell's equations and the elasticity equations. For Maxwell's equations the null space 
of C contains all symmetric matrices, while for elasticity the null space of C contains 
antisymmetric matrices, which is a space of lower dimension. This manifests itself in the 
different boundary conditions: at an interface u is required to be continuous, while only 
the tangential component of E is required to be continuous. In this respect Maxwell's 
equations have some similarity with the acoustic equations, which (see (12.51) ) take also a 
form analogous to (II. 2p or (11.41) . with 

and the null space of C contains all matrices which have zero trace, and only the normal 
component of u is required to be continuous at an interface. 

There are also linguistic differences when one discusses elasticity compared to elec- 
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tromagnetism. When one wants to study equation (1.2) in a bounded domain, boundary 
conditions are needed. A natural condition is to fix the value of n ■ C(9E/(9x where n is 
the external unit normal to the boundary of the domain. The analog boundary condition 
in the elastodynamic case is well known and called the surface force F applied to the 
medium. Thus we will call the applied surf ace free current the value of (2Ci;)~^n- C(9E/9x 
on the boundary and denote it J. (The additional factor of {iuj)~^ is introduced because 
ia;j in (11.21) plays the role of f in (11.41) .) This is not a usual way of speaking in the electro- 
magnetic framework as the value of {iu!)^^n - C(9E/(9x is nothing else than the tangential 
part of the magnetic field H at the boundary. The interest of such a vocabulary appears 
later. 

When one wants to study equation (1.2) in a domain Q which is divided in two sub- 
domains Qi, Q2, one has to write jump conditions on the dividing surface. This condition 
is the continuity of the tangential part of H (the analog of which is the continuity of the 
normal part of the stress in elastodynamic framework). Alternatively one says in the 
elastodynamics framework that fl2 exerts on Qi a surface force F while Qi exerts on Q2 
the opposite surface force — F. This action-reaction law makes the link with the separate 
study of both subdomains as it fixes the needed boundary conditions for these studies. 

In a similar way we can say that Q2 exerts on Qi a surface free current J while 
fli exerts on Q2 the opposite surface free current —J. This way of thinking needs some 
practice to become natural and the reader should be aware that this formulation does not 
mean, in any way, that there exist actual free currents in the material (just like action- 
reaction law does not imply the existence of actual surface forces inside the domain). 
However the surface force F that fl2 exerts on Qi has an equivalent effect on Qi as a 
body force f concentrated at the boundary replacing the stress field in Q2, and similarly 
the applied surface free current J that ^2 exerts on Qi has an equivalent effect on Qi as 
a free current j concentrated at the boundary replacing the H field in 0,2- 

It is well known that when j = the Maxwell system of equations remains unchanged 
when one interchanges the roles of E and H and of e and fi. Therefore for each elec- 
tromagnetic circuit discussed here, there is a dual magnetoelectric circuit (ME-circuit) 
obtained by making these replacements. Instead of speaking about applied free surface 
electric currents, we could speak about applied free surface magnetic monopole currents. 
These are then truly unphysical, but their introduction is again just a device for keeping 
track of boundary conditions. 

Throughout the paper we use infinite or zero values of various moduli. From a physical 
viewpoint one should think of such moduli as just being positive and real and extremely 
large or extremely small. From a mathematical viewpoint one should think of taking the 
limit as these moduli approach infinity or approach zero. Generally values of the permit- 
tivity and permeability near zero or infinity are difficult to achieve. However, using res- 
onance effects (ISchelkunoff and Friis 19521 Pendry, Holden, Robbins, and Stewart 1999 ) 
very small or very large values which are almost real and positive may be achieved over a 



narrow frequency range. The importance of this was recognized by Engheta, Salandrino, and Alu (2005) 
and Engheta (2007) who realized one could build nanoscale equivalents of electrical cir- 
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cuits using such materials: a material with e near zero electrically insulates the circuit, 
while a material with e near infinity provides the necessary electrical connections. 

From a physical viewpoint perhaps the greatest barrier to the construction of electro- 
magnetic circuits is the use of a matrix, which is the electrodynamic equivalent of a void 
in elasticity, and has an extremely large value of the magnetic permeability /i, and (al- 
though it is not clear it is necessary) an extremely small value of the electric permittivity 
e. (In the case of the dual circuits, one would need the reverse). In fact it is not necessary 
that the circuit be embedded in a body with these properties, only that a material with 
these properties clads the circuit. Also note that the Maxwell equations (11. 2p remain 
valid if one divides /^(x) everywhere by a constant k and correspondingly multiplies j(x) 
and £(x) by k. Therefore it should be possible to renormalize the moduli in the EM- 
circuit in such a way, that the moduli in the matrix take more realistic values, perhaps 
even that of empty space with e = ^ = 1. This is similar to the way a spring- mass 
network can still function when embedded in an elastic material provided the springs are 
appropriately stiff, the forces sufficiently strong, and the masses are sufficiently heavy. 

The objective of this paper is to introduce the concept of EM-circuits and their basic 
properties. The approach is formal, but will hopefully motivate future analytical and 
numerical work to place the treatment given here on a firm foundation. 



2 Transverse electric EM-circuits 

We are interested in the Helmholtz equation 

V ■ (1/yuVE) = -euj^E, (2.1) 

describing three-dimensional TE electromagnetic wave propagation, where uj is the (fixed) 
frequency of oscillation, E(x) = (0, 0, E) is the electric field, £:(x) is the electrical per- 
mittivity, /i(x) is the magnetic permeabihty, and all of these quantities do not depend 
on x^. Given the electric field component E the associated in-plane magnetic field is 

H = -[V(/Ucj)]R^VE, (2.2) 

where 

Rx^C; (2.3) 

is the matrix for an anticlockwise rotation by 90° in the plane. When in some subregion 
f2, the moduli are real, positive, and do not depend on frequency, the electromagnetic 
energy stored in Vt is 

W{n) = [ [fi\li\' + e\E\']/A. (2.4) 
Jn 

We are only interested in solutions such that Vr(fi) remains bounded in all subregions Q 
where the moduli tend to zero or infinity, remaining real and positive in this process. 
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The Helmholtz equation (12.11) is mathematically analogous to the acoustic equation 
which in two-dimensions reads as 

V ■ (1/pVP) = -{1/k)u:^P, (2.5) 

where -P(x) is the pressure, fi;(x) is the bulk modulus, and p(x) is the density. Given the 
pressure field -P(x) the associated velocity field of the fiuid is 

v = -[V(p^)]VP, (2.6) 

and when in some subregion VL the moduli are real, positive, and do not depend on 
frequency the time averaged elastokinetic energy in VL is 

W{n)= [ [p|v|2 + |P|V«:]/4. (2.7) 
Jn 

(The extra factor of 2 arises because the physical velocity and pressure is the real part 
of e-*^*v and C'^^P). 

To conceive TE electromagnetic circuits we just have to understand how discrete 
acoustic networks are made and transcribe their structure in terms of electrodynamic 
quantities. As illustrated in figure [2] we consider a network of channels connected by 
junctions. Each channel has parallel sides and contains a segment of incompressible, 
non-viscous, fiuid with some constant density p possibly varying from channel to channel, 
moving in a time harmonic oscillatory manner in response to time harmonic pressures at 
the junctions. We define the entire cavity associated with a junction to be the cavity at 
the junction, plus the remaining region in the channels not occupied by the incompressible 
fluid. If the junction is a terminal so that there is an open channel leading to it, we insert 
a segment of massless incompressible fluid in that channel to keep track of the response 
of the acoustic network. 

Each entire cavity contains a compressible, non-viscous, massless fluid with com- 
pressibility possibly varying from cavity to cavity. In this model the compressibility is 
localized in the cavities and the mass is localized in the channels between cavities. The 
surfaces between the compressible and incompressible fluids are assumed to be flat and 
perpendicular to the channel. When one conceives an acoustic network, the area outside 
the network is assumed to be occupied by a rigid boby (or cladded by rigid tubes). How- 
ever we can, in a equivalent way, assume that this area VLq is fllled by an incompressible 
fluid having inflnite density, i.e. with k = p = oo. Indeed, the inflnite density, and the 
boundedness of W{flo) ensures that the velocity v (as one might physically expect) will 
be zero outside the network and consequently that the walls of the network remain flxed. 
The incompressibility ensures that both sides of (12.51) are zero in the matrix, without re- 
quiring that P = in the matrix. Hence the acoustic equation (12. 5p desribes the system 
in the whole space. 

In each entire cavity, Qj, where 1/p is inflnite (or more precisely is real positive and 
approaches inflnity), VP must be zero, since otherwise p|vp = |VPp/(pcj^) would be 
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Figure 2: A two-terminal discrete acoustic network. In the idealized model the four 
cavities contain compressible massless fluid, while the grey shaded fluid plugs in the five 
tubes contain incompessible fluid with some mass. 

infinite and W{Qj) would be unbounded. Thus, as expected, the pressure is constant in 
each junction region. Within each segment of incompressible fluid of constant density, 
both sides of (12.51) must vanish, which implies V^P = in each such segment. From 
the boundary conditions (that P is constant at the ends of the fluid segment, and at the 
sides n • VP = 0, since n ■ v = 0) it follows that the pressure P will be constant in each 
cross section normal to the channel, and will vary linearly along the fluid segment. The 
fluid velocity v will therefore be constant in the segment, and directed parallel to the 
fluid channel. From (12.61) we see that in a channel joining cavity i and cavity j the fluid 
velocity in the direction of the channel, from j to k will be 



where ijk is the length of the fluid segment, pjk is its density and Pj and Pk are the 
complex pressures at junctions j and k respectively. This is basically Newton's second 
law, relating the acceleration of the fluid segment, —iuvij, to its mass and the force 
acting on it. 

In entire cavity j, (12.51) and (12. 6p imply 



which when integrated over the entire cavity implies, by the divergence theorem. 



Vjk = -[i/{pjk^^ijk)]{Pk - Pj), 



(2.8) 



V- V 



{i/ k)ujPj 



(2.9) 



m 




(2.10) 



k=l 
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where aj is the area of the entire cavity, Ki is the bulk modulus of the fluid within it, and 
we have assumed that m channels enter the cavity, each with width h and carrying a fluid 
segment with velocity Vjk- This is essentially Hooke's law, applied to the compressible 
fluid occupying the entire cavity. If the cavity is a terminal and there is an open channel 
carrying a current Ij into it, then the relation fl2.10p takes the modified form 

h 

Ij + ^ hvjk = ajiil Kj)ujPj. (2.11) 

k=l 

Given the pressures Pj at the terminal cavities, the equations (12.80 and (12.101) provide a 
discrete set of equations, which can be solved for the pressures in the other cavities, the 
velocities of the fluid plugs in the channels, and the currents Ij flowing into the terminal 
cavities. Assuming that the cavities are numbered in such a way that the first n are 
terminals, and the remaining ones are not, the response of the network is expressed in 
terms of the linear relation 

I = MP (2.12) 

between the set of pressures P = {Pi, P2, ■ . . , Pn} at the terminals and the currents 
I = {/i, J2, . . . , /„} flowing into them. 



Everything carries through to the electromagnetic case where the fields are transverse 
electric (TE). By comparing (12. 5p and (12. 6p with (12. ip and (12. 2p we see that fi and e 
play the role of p and l/zt; E plays the role of P and H plays the role of R_lv (and 
is therefore perpendicular to the channel walls). The channels themselves are now thin 
plates containing a material with e = and 7^ 0. The cavities are now aligned dielectric 
cylinders with e ^ and fi = 0. The electric field is constant in each cylinder, which 



also can be seen directly from the result of Silveirinha and Engheta (2006) who show 
that, for the dual transverse magnetic (TM) problem, cylinders having 7^ and e = 
have a constant magnetic field in them. It follows that H = in the matrix by direct 
analogy with the acoustic case where v = in the matrix. (Although magnetic fields 
tend to be concentrated inside a material with positive and very large permeability, this 
concentration refers to the B field and not to the H field). 

We call such a circuit a transverse electric EM-circuit (see figure [3]). Each equation 
we discussed has its analog. For example (12.111) becomes 

m 

Ij + hHjk = ajiujEjEj, (2-13) 

k=l 

where Hjk is the value of H in the direction perpendicular to the walls of the plate jk, Ij 
is the line integral of H across the open channel, Ej is the dielectric constant of cylinder 
j while Ej is the electric field in the cylinder j. So the left hand side of (12.131) is the line 
integral of H around the terminal dielectric cylinder, while the right hand side of (I2.13P 
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Figure 3: A two-terminal transverse electric EM-circuit which is the exact analog of the 
discrete acoustic network of Figure 1 

is the total displacement current flowing through the cylinder. Thus f l2.13p is nothing 
but Ampere's circuital law (with Maxwell's correction). Notice that instead of having an 
open channel one could have a free current —Ij flowing next to the terminal dielectric 
cylinder. 

The analog of ( 12. Sp is 

Hjk = -\l/{^ljM,k)]{Ek ~Ej), (2.14) 

where iXjk is the permeability of the plate jk. This is Faraday's law of induction, relating 
the time derivative of flux of B through any rectangle with two opposite sides along the 
dielectric cylinders j and /c, to the line integral of E around this rectangle. Since H is 
constant and perpendicular to the plate walls, it follows from (12. 2p that E in the plate 
depends linearly on xi and X2 in such a way that it is constant perpendicular to the plate. 

In the particular case when the cylinder j has zero dielectric constant, i.e. Ej = 0, (so 
that the junction is the analog of a cavity filled with incompressible fluid) (12.130 becomes 

m 

+ 5^ = 0, (2.15) 

k=l 

where Ij = if the cylinder is not a terminal cylinder. If all cylinders have zero dielectric 
constant, then we call the circuit a transverse electric magnetic circuit (M-circuit). 

It is now important to understand how can these circuits be used and in particular 
how they can interact with ordinary materials. The problem is analogous for connecting 
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an acoustic discrete network to an ordinary acoustic three-dimensional domain. Assume 
that the matrix with fi = oo and 6 = only has finite extent, and is surrounded by 
space with /i = e = 1, in which there are TE fields. Also suppose each terminal edge 
is connected to the exterior by an open channel, of width h, containing material with 
e = fi = 0. The external field E (which is the anolog of the pressure) will fix the mean 
value of Ej at every open channel mouth j. Then the response of the transverse EM 
circuit will determine the values Ij, that is of H (which is analog to the velocity) at 
the open channel mouths. Let Qq denote the region occupied by the circuit plus the 
remaining matrix. On the part of the boundary of Qq which corresponds to the matrix 
we have H = 0. Hence the ii^-to-tangential value of H map (which is equivalent to the 
Dirichlet to Neumann map) of Qo will be governed by the response of the circuit. It will 
be completely different from a pure matrix (for which the tangential value of H vanishes 
on dQo) or from void where e = fi = 1. 

Note that the external field E will fix the value of the electrical field at each open 
channel mouth in a efficient way if h is large enough. If h is too small this connection 
will be weak, and E near each mouth will be strongly affected by flux of R^H (which is 
analogous to current in the acoustic setting) through the narrow channel openings. (In a 
region near but not too close to each mouth the E field will be like that generated from 
a line source.) However this problem can be corrected by adding at each open channel 
mouth a material with e = fi = (see figure H]). Hence the value of E at each mouth will 
be fixed by the value of E on the 'relatively large' cap boundary, and the fiux of R^H 
through the channel will be transferred to the outside of the cap. This is similar to the 
way Silveirinha and Engheta (2006) has, for the dual TM problem, suggested the use of 
materials with e = for transfering energy through narrow openings. 



3 Electromagnetic circuits in the general case 

We need to generalize the EM-circuits to allow for fields that are not transverse electric. 
Like in the transverse electric case the circuit will be composed of plates of material 
having e = and /i 7^ 0, joined by cylinders of dielectric material having e 7^ and 
/i = 0, embedded in a matrix having fi = 00 and e = 0, so that the matrix is the 
electrodynamic equivalent of a void in elasticity according to (ll.2p - fll.4p . 

The plates and the dielectric cylinders play the physical role in our circuits that 
springs and masses play in an elastic network, despite the fact that they are completely 
different geometrical objects. The assumption that the electromagnetic energy density in 
the matrix remains bounded as /i 00 and e ^ when the moduli fi and e are positive 
and independent of frequency, again implies that H = D = in the matrix. (Note that 
if H = in the matrix then necessarily D = zV x H/cj is also zero). 

We emphasize that when e and fi are real and positive in the matrix and fi is very 
large, while e is very small then there certainly exist (high energy) solutions where H 
in the matrix is not small: after all an electromagnetic wave could propagate there, and 
its amplitude scaled as one desires. However, we believe (and this needs to be rigorously 
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Figure 4: A semicircular plug of material with e = /i = can serve to couple the open 
channels of an transverse electric EM-circuit with exterior TE fields, allowing the H field 
to be transfered to the exterior with negligable drop in E. 

verified) that the solutions in the matrix almost decouple from the solutions in the elec- 
tromagnetic circuit when /i is very large and e is very small. This should be similar to 
the way electromagnetic fields almost decouple at a planar interface between two non- 
absorbing media, 1 and 2, for which there is a large mismatch in the electromagnetic 
impedances rji = fii/ei and 772 = when 771/772 is very large then a plane electro- 

magnetic wave incident from either side of the interface will have only a tiny portion of 
its energy transmitted. 

Alternatively, and as kindly suggested to us by a referee, one may assume that in the 
matrix the product e/i is large and negative. Then electromagnetic fields in the matrix 
will be confined within a small skin depth of the surface which tends to zero as e/i —00, 
again implying that H = D = in the limit as /i — > 00 and £ — in such a way that 
e/i — s> —00. 

Let us now analyze in detail the response of each plate. The plate could be polygonal 
in shape, but for simplicity we use a basic element which is a very thin triangular prism, 
of uniform height h containing a material having e = and /i 7^ 0, the top and bottom 
faces of which are surrounded by the matrix, as illustrated in figure [51 We call this 
element a magnetic element. 

Let us choose our coordinate system so the bottom surface of the prism is at X3 = 
and the top surface at X3 = h. The triangle at the bottom of the prism has vertices 
q, r and s, labeled in an anticlockwise order when viewed from above, and edges qr, rs 
and sq. Let fiqrs denote the constant value of fj, within the prism. Since V ■ H = and 
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X 



Figure 5: A triangular prism containing material with e = and jj, is one basic 
element of a magnetic or electromagnetic circuit. Here the arrows denote direction of the 
magnetic field H within the prism. 

V X H = in the prism it follows that H = Vtp where Atp — 0. Also since the tangential 
component of H = V^j is zero at the top and bottom surfaces of the prism, it follows 
that ip is constant on the top and bottom plates: ip is like the potential between two 
closely spaced capacitor plates. Hence H is essentially constant within the prism and 
normal to the top and bottom surfaces, i.e. H = (0, 0, H^) where cannot depend on 
X3 since V • H = within the prism. (In fact H will only be approximately constant due 
to fringing fields which, however, should become negligible away from the edges, in the 
limit as the prism becomes very thin.) 

Let x' = {xi,X2) denote coordinates in the plane. Assuming the point r is at xi = 
X2 = 0, the three edges of the triangle lie along the three lines 



each parameterized by a where is the unit vector directed from vertex m to vertex n, 
and Xq is a point along the edge sq. In electromagnetic circuits we constrain the tangential 
component of the electric field to take constant values Egrtqr, Erstrs, and Egqtgq along 
the three sides gr, rs and sq of the triangle. (As we will see later the presence of dielectric 
cylinders along these edges will allow this constraint to be satisfied). Let l^s and Igq 
denote the lengths of the edges qr, rs and sq. Then Faraday's law of induction applied 
to a circuit around the triangle implies 



x' = at, 



x = at, 



x' = Xq + at, 



(3.1) 




(3.2) 
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where c = iujjJiqrsH^ and a^^s is the area of the triangle. 

To find an exphcit expression for the electric field in the prism, although it is not 
clear we need it, let us assume that e in the prism is arbitrarily small, but non-zero (and 
many factors greater than the e in the matrix which we treat as being zero), so that in 
the prism V ■ D = implies V ■ E = 0. Since D = in the matrix and V ■ D = it 
follows that D and hence E are tangential to the top and bottom surfaces of the prism. 

Having a material with zero permittivity outside the prism allows us to have a non- 
zero E field there. Then the tangential components of E can be continuous across the 
top and bottom surfaces of the prism. It is not clear that this zero permittivity in the 
matrix is necessary. One could instead have e 7^ and E = outside the prism, with a 
concentrated surface B current to compensate for the jump in the tangential component 
of E across the surface. Such a concentrated surface B current should be allowed since 
/i = 00 in the matrix. (Similarly in an elastic network, it is not necessary that the 
surrounding material have density p = 0, although that is the case when the surrounding 
material is void. If p is non-zero and C is close to zero then only a small boundary layer 
near the elastic network will move.) 

Since V x E = zcu/xH, we infer that 

E = (-cx2/2,cxi/2,0) + V0, (3.3) 

where A0 = and without loss of generality one can assume that = at the origin. 
The potential (f) satisfies the Neuman boundary conditions that n • V0 = on the top 
and bottom surfaces of the prism, and Dirichlet boundary conditions on the sides of the 
prism (specifying the tangential value of V0 around the sides, and the value = at 
the origin determines the value of along the sides). Thus is uniquely determined and 
a simple calculation using (13. ip . the identity 

^qr^qr ~t~ ^rs'^rs ~t~ ^sq^sq ('^*^) 

(as follows from the fact that the edges form a triangle) and the fact that agrs = 
isqtsq ■ R-±xo/2 is the area of the triangle (as can be easily seen by choosing Xq to be per- 
pendicular to tgq ) shows that the boundary conditions are satisfied with V0 = (ai, 02, 0) 
where a = (ai, 02) is constant and determined by 

tqr ■ a = Eqr, ^rs ' ^ ~ Ej-g. (3-5) 

(The condition (13.21) ensures that E ■ t^g = Egq along the edge sq). 
It is natural to introduce three new variables 

Vqr Eq^Eq^^ ^rs ^rsErs^ ^sq ^sqEgq^ ('^*^) 

which when i/3 = would represent the potential drops along the three edges. Then 
(13.21) implies that Hy, depends on Vqr, K-s, and Vgq only through the sum Vqr + + Vsq- 
In keeping with the vocabulary introduced in the introduction, the material surround- 
ing the edges of the basic element exerts total applied surface free currents J^^, J^^ and 
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Jlq along the edges qr, rs and sq, flowing in directions tqr, trs, and t^g, where the su- 
perscript is kept to signify that the currents are associated with the triangle qrs. (Here 
total signifies that these are the applied surface free currents integrated over the width 
of each edge, but from now on this will be assumed so we will drop the word total). In 
other words, the boundary conditions on the edges of the basic element are essentially 
the same as if we completely surrounded the basic element by matrix material with e = 
and /X 7^ having H = and inserted these surface free currents along the edges. 

These currents are all equal, and from Ampere's circuital law applied to a circuit 
around each edge take the value hH^, by virtue of the fact that H is constant within the 
triangular prism. It may seem superfluous to keep track of the three currents J^^, J^^ and 
JJ^ when they are all equal. However, consider the analogous elastodynamic framework: 
to write the balance of forces at each node, one introduces the forces that each spring 
exerts on each node even though the forces exerted by a spring on its two extremity 
nodes are equal and opposite. Without introducing J^^, J^^ and JJ^ it would be difficult 
to derive an expression for the response matrix of an general electromagnetic circuit, as 
we do in section 5. 



Thus we have the relation 





(3.7) 



where kqrs = h/ {aqrsl^qrs)- We use the quantities T rather than the free currents J so that 
the matrix entering the above relation is real and so that the parallel with elastodynamics 
is maintained, since ituj in (11. 2p plays the role of the body force f in (II. 4p . 

Also our introduction of the variables V rather than the variables E ensures that the 
matrix is symmetric which is desirable since this property will then extend to the matrix 
describing the response of electromagnetic circuits with many elements. (Again, this is 
why it is important to introduce the three currents J^^, J^^ and JJ^ rather than just a 
single current.) This relation (13.71) . which is essentially Faraday's law of induction, is the 
analog of Hooke's law. 



^ / ^rs '<y ^rs ^rs '<y ^rs j / j 



describing the response of a spring, where and F^ = — F^ are the forces node r and 
node s, respectively, exert on the spring joining these nodes, [which is the opposite of 



the definition given in Milton and Seppecher (2008)], and are the displacements 



at these two nodes, trs is the unit vector pointing from node r to node s, and krs is 
the spring constant. Note that the matrix entering both relations (13.71) and (13. 8p is real, 
symmetric, degenerate and positive semidefinite. Also iuJ is playing the role of a force, 
and 1/fiqrs is playing the role of the elastic spring constant (to within a proportionality 
factor) as might be expected by comparing (II. 2p and (II. 4|) . 
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Figure 6: A magnetic circuit (M-circuit) is obtained by joining together a collection 
of triangular prisms, of the type illustrated in figure [5l At each edge there is a small 
diameter cylinder, not illustrated, having e = fi = 0. Here the terminal edges are marked 
by thicker lines. An EM-circuit is obtained when a selection of the cylinders along the 
edges are assigned are a dielectric constant e 0. 
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A magnetic circuit (the analog of a elastic network with springs but no masses), as 
illustrated in figure El is a collection of such triangular prisms, joined at common edges 
by cylinders having n = e = and with a constant diameter d of the order of h. In fact it 
is desirable to take e in these cylinders arbitrarily small but non zero, since then E will be 
constant along the cylinder because D is (essentially) constant. Edges in such a magnetic 
circuit (M-circuit) play the role of nodes in an elastic network, and just as applied forces 
are confined to the terminal nodes in an elastic network, so too can applied free currents 
be confined to a subset of the edges in a magnetic circuit. We call these the terminal 
edges, and we call the others internal edges. If a magnetic circuit contains an internal 
edge qr which is connected to only one triangle qrs, then T^^ = T^^ = TJ^ = and the 
triangle qrs can be removed without effecting the response of the network. (Analogously, 
if a spring network contains an internal node with only one spring and no mass attached 
to it then that spring can be removed without affecting the response of the network). 
Thus we can restrict attention to magnetic circuits where all internal edges are connected 
to at least two triangles. 

Consider an internal edge qr where m triangles meet at a cylinder. Since T^^/ (iuh) is 
the value of the constant magnetic field H within the triangular prism qrs, at an internal 
edge qr where m triangles meet at a cylinder, we have 

m 

^r;, = 0, (3.9) 

as follows from Ampere's circuital law that the line integral of H around the cylinder is 
zero. Equation (13.91) is analogous to the balance of forces at a node in a spring network: 
the sum of all free currents must be zero if there is no net free current. At a terminal 
edge qr, Ampere's circuital law implies 

m 

J2T^r = Aqr, (3.10) 
s=l 

where Agr/ (iuj) is the free current applied to that edge. 

We label the edges in the network so that no edge is repeated twice, i.e. if qr labels 
an edge in our list, then the label rq does not appear in the list. This essentially assigns 
an arrow (from q to r) to each edge, and it may be impossible to assign arrows so that 
no two arrows point to the same vertex in every triangle in the circuit. Accordingly, for 
example, we may want the relation (13. 7p to involve T^^. and Vsr rather than T^^ and Vrs 
when the label rs does not occur in the list. To eliminate the unwanted variables in (13. 7p 
we can then use the relations 

Tr's = -T!r, Vrs = -V,r, (3.11) 

which hold for all r, s and q. Thus (13. 7p becomes 

Tl = ILUJ^, I = kgrs 1-1 1 -1 I I I , (3.12) 



1 
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Figure 7: An electromagnetic circuit (EM-circuit) has small diameter dielectric cylinders 
with £ 7^ and fi = along a selection of edges including possibly the terminal edges. 
Here there are three terminal edges are marked by thicker lines. The dielectric cylinders 
should be much thinner than drawn here. One internal edge, marked by the thin line, 
has a cylinder with e = fi = attached to it. 

and still involves a real, symmetric, degenerate, positive semi-definite matrix. 

To obtain an electromagnetic circuit from a magnetic circuit (including those magnetic 
circuits where some internal edges are only connected to one triangle) we assign a non- 
zero value to the dielectric constant (to some or all) of the cylinders, of diameter d, at 
the junctions of the triangles. (This is analogous to adding mass to the nodes of a spring 
network) An example is illustrated in figure [71 At any vertex where two or more dielectric 
cylinders meet we need to make sure there is a good electrical connection between the 
dielectric cylinders to allow displacement current to fiow between the cylinders. 

Now at an internal edge qr where m triangles meet at a dielectric cylinder the junction 
locally looks similar to the junction in a transverse electric EM circuit where m plates 
meet at a dielectric cylinder, and so one expects an equation similar to 02.131) to hold. 
Ampere's circuital law (with Maxwell's correction) taken around a circuit surrounding 
the cylinder qr implies 

m 

Y,K=^'9,rVgr, (3.13) 

s=l 

where gqr = ncPegr/ (4£gr), in which Sgr is the dielectric constant of the cylinder. The term 
on the right arises from the fact that —i{'KcP/A)ujeqrEqr is the total displacement current 
flowing through the dielectric cylinder. Inside the cylinder D = Vx (since V ■ D = and 
V X E = 0) where n • Vx = at the cylinder walls (since D = in the matrix and in the 



18 



triangular elements). At the cylinder ends one has some flux of D. From the solution to 
this Neumann problem D will be essentially constant inside the small diameter cylinder 
away from the ends. This justifies our assumption that the electric field takes constant 
values along the edges of a magnetic triangular element, at least when there are dielectric 
cylinders along each of these edges. 

The equation (13.131) which is essentially the same as (12.131) when Ij = 0, is the analog 
in an elastic network of Newton's law, 

m 

J2F'^ = U^TJlrUr, (3.14) 
s=l 

describing the motion of a mass rrir at a node r where m springs meet. At a terminal 
edge (I3.13P generalizes to 

m 

Y,T-^r = Aqr+UJ^9,rV,r (3.15) 
s=l 

where again Aqr/ iiuj) is the free current applied to that edge. The equations (13. 7p . (13.131) . 
and (13.151) hold for each triangle and each edge, and provide a system of equations which 
can be solved to determine the response of an arbitrary electromagnetic circuit. This will 
be done in section 5. 

The mathematical idealization of an electromagnetic circuit is obtained by taking 
the limits /i ^ and ^ 0, while say keeping the ratio d/h fixed. The moduli of the 
constituent materials need to be scaled in such a way that the parameters entering the 
final equations, such as kq^s and Qqr, remain fixed. Thus one should take /Xg^s proportional 
to h (and thus very small) and Sqr proportional to l/c?^ (and thus very large). 

4 Acting upon an electromagnetic circuit and creat- 
ing virtual free currents 

One might ask how one could conceivably act on an EM-circuit, and measure its response. 
A possible scenario, as sketched in figure [HI might be to be to have electromagnetic fields 
incident on a body, say a cube, of material with n = oo and 6 = containing an EM- 
circuit, with no two terminal edges sharing a common vertex, positioned in such a way 
that only the terminal edges are exposed at the surface of the cube. Let us suppose that 
there are no dielectric cylinders attached to the terminal edges. Then the electric field 
will not be constant along each terminal edge. If qr is a terminal edge between points q 
and r both on the same face of the cube, then Faraday's law of induction implies that 
the line integral of E along that terminal edge will play the role of the quantity Vqr in 
the electromagnetic circuit so that (13. 7p remains satisfied. 

The EM-circuit causes the magnetic field H outside the body to be altered in such a 
way that Ampere's circuital law (with Maxwell's corrections) holds around each terminal 
edge. If one was not aware of the existence of the EM-circuit, from outside the body 
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|lI=1, 8=1 



|Ll=oo, 8=0 



Figure 8: A cube of material with n = oo and e = containing an EM-circuit with four 
terminal edges exposed on the faces of the cube. From the outside it will look as if there 
are free currents flowing along the terminal edges (although in reality they do not exist) 
with the endpoints acting as sources and sinks for the displacement current field outside 
the cube. These virtual free currents will be generated according to the values of the 
line integrals of E along the terminal edges, and according to the response matrix of the 
EM-circuit. In the dual setting, it will look like a ME-circuit generates virtual magnetic 
monopole currents. 
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it would look as if the magnetic field H near the terminal edge qr was generated by a 
free-current —Aq^/ {iuj) flowing from q to r. In other words, if one incorrectly assumes 
that H = throughout the cube, then Ampere's circuital law would falsely imply the 
existence of this free current, which we call a virtual free current, flowing along the 
terminal edge. It is nothing else but the surface free current which the EM-circuit exerts 
on the surrounding material at the terminal edge qr. (In a similar fashion one can insert 
a mass-spring network into a cavity in an elastic body, with only the terminal nodes 
attached to the boundary of the cavity. If one is not aware of the existence of the spring 
mass network from outside the cavity, it would look like the stress field in the body was 
altered by concentrated forces acting at the positions of the terminal nodes.) 

Now the internal edges will carry some displacement current Aqr/{ii^) out of the 
vicinity of the point q and a displacement current Aqr/ {iui) into the vicinity of r. If one 
is not aware of the existence of the electromagnetic circuit it would look like the point p 
is a current source and the point r is a current sink: it would look like the ends of the 
virtual free-current ~Aqr/{iuj) along the terminal edge, act as sources and sinks for the 
displacement current outside the body. 

If the thickness h of each terminal edge is very small, then the coupling between 
the electromagnetic circuit and the fields in the exterior will be weak. As in the case 
of transverse electric EM-circuits, small virtual free-currents along the terminal edges 
will cause the field E(x) to be modified in the near vicinity of each terminal edge. One 
suggestion to enhance the coupling is cap each terminal edge qr with a e = fi = 
semicircular cylinder of length £qr and diameter do, where do is not small. At the two 
ends of this cylinder one could attach e — oo, fj, — quarter spheres of diameter do, to 
allow the displacement current to enter and exit the points q and r with little resistance. 
In these quarter spheres E = 0. Faraday's law of induction then implies that the line 
integral of E along the outer surface of the semicircular cylinder will equal the line integral 
of E along the terminal edge. 

5 A formula for the response matrix of an EM-circuit, 
and the properties of this response matrix 

In a magnetic or electromagnetic circuit with n terminal edges let us suppose these edges 
have been numbered from 1 to n. Then the response of the network is governed by the 
linear relation 

A = WV (5.1) 

between the terminal variables A = {Ai, A2, . . . , An) which measure the real or virtual 
free currents at these edges, and the variables V = {Vi, V2, . . . , V,i) which measure the 
line integral of the electric field along these edges. When all the edges in the circuit are 
terminal edges the response matrix W equals a symmetric matrix W° with an especially 
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simple form. From (13.71) and f l3.15p the diagonal elements of W° are given by 

W^r,,r = -^^9<ir + J2Krs, (5.2) 

s 

where the sum is over vertices s such that qrs is a triangle in the circuit, while the off- 
diagonal elements W^g^st ^^^^ when qr and st are not two edges of some triangle in 
the circuit, and the remaining off diagonal elements are each given by one of the formulas 

TI/O _ Tj/0 _ L Ti/0 _ Ti/0 _ _L /e- o\ 

'^'^qr,rs qr,sq '^qrsi qr,ST qr,qs '^qrsi 

according to what edge labels are in our list, where qrs is a triangle in our circuit. 
Suppose we divide these edges into two groups, and order the edges so that one group 
comes first. Then the matrix relation fl5.ll) takes the block form 



(5.4) 



where Ai and Vi are the set of variables associated with the first group and A2 and V2 
are the set of variables associated with the second group. Now consider the case where the 
first group are terminal edges, while the second group are internal edges. Then A2 = 
and (15. 4p implies Ai = WVi with the response matrix W of the circuit being the Schur 
complement 

W = W?, - W?2(W°2)-^(W?2)^. (5.5) 

This is our formula for the response matrix of an arbitrary electromagnetic circuit. In 
particular it shows that the response matrix is always symmetric. It may be that the 
matrix W22 is singular, in which case if Ai is finite there are generally restrictions on 
the possible values that Vi can take. 

Also recall that the matrix entering the relation (13.71) is positive semidefinite. There- 
fore if ^qrs has a non- negative imaginary part, and hence kqj-s has a non-positive imagi- 
nary part, for each triangle qrs in the circuit and e^^, and hence gq^.., have a non-negative 
imaginary part for each edge qr in the circuit, the imaginary part of W° will be negative 
semidefinite, being a sum of negative semidefinite matrices. It follows that the quantity 

S = A!-Y"- A" ■ V = -V ■ (W°)"V' - V" ■ (W°)"V" (5.6) 

will be non-negative, where the primes denote real parts, and the double primes imaginary 
parts. In particular if A2 in (15.41) is zero, the left hand side of the above equation reduces 
to 

s = A[- Y'l - A'[ ■ v; = -v; ■ w"v; - v;' ■ w"v;', (5.7) 

and since this is non-negative for all values of Vi we deduce that W", like (W°)", is 
negative semidefinite. 
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6 The energy stored and dissipated in an electro- 
magnetic circuit 

To obtain a formula for the time averaged energy stored in an electromagnetic circuit 
let us assume the moduli in the electromagnetic circuit are real and do not depend on 
frequency. Recall that the physical electric and magnetic fields are the real part of Ee~*'^* 
and He~*^*. Locally the time averaged electric and magnetic field energy densities will 
therefore be £|Ep/4 and /i|Hp/4. In the magnetic element qrs discussed at the beginning 
of section 3 the time averaged stored magnetic energy will be 

ha,rs^l,rs\H^\V^ = \T',r?/{^^\rs). (6.1) 

while in the dielectric cylinder qr the time averaged stored electrical energy will be 

{T^l,r(f/^)e,r\E,r?/^ = |K,r|'v/4- (6-2) 

The total electromagnetic energy stored in the electromagnetic circuit will be a sum of 
such expressions taken over all magnetic elements and dielectric cylinders in the circuit. 
An appealing feature is that the resultant expression only depends on uo and the parame- 
ters T, V , k and g characterizing the electromagnetic circuit, and not on the parameters 

Now let us consider how much electromagnetic energy is dissipated into heat within 
the electromagnetic circuit when the moduli are complex and depend on frequency. Lo- 
cally the time averaged electrical and magnetic power dissipated into heat per unit volume 
will be ue"\Fj\'^/2 and a;/i"|Hp/2, respectively. Within the magnetic element qrs this will 
integrate to 

u;H,,/i;'-l^3lV2 = il/k,rs)"\T^r\Vi'^^) 

= mr)\TyKrsr - iT:r)"K-/Krs)Vi2uj). (6.3) 

Now we can substitute (13.71) into this, and associate a portion of the resultant expression 
to each edge, where the portion assigned to edge qr is 

mr)K - (^,;)X]/(2^)- (6-4) 
In the dielectric cylinder qr the time averaged electrical power dissipated into heat is 

{nigrdVA)uje';^\E,r\V2 = uj\V,r?glJ2 

= [{-u^gMXr - {-^^gM'Xr]K'2^)- (6-5) 

Adding up all the contributions (16. 4p and (16. 5p associated with edge qr and using the 
relation (I3.15P we see that the total contribution associated with edge qr is zero for an 
internal edge and 
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for a terminal edge. By summing this expression over all terminal edges we see that the 
quantity S/{2u), where S is given by fl5.7l) . is the time averaged electromagnetic energy 
converted into heat in the circuit. 

This is consistent with Poynting's theorem. Suppose we attached to the edge qr a 
rectangular plate of thickness h and width igj. in which there is a magnetic field H with 
component Hqr = —Agr/ {iuh) perpendicular to the plate (and 

surrounded by material with n = oo and e = 0) so that Ampere's circuital law (with 
Maxwell's corrections) is satisfied around the terminal edge. At any instant in time the 
flux of energy into the terminal qr will be hiqr{Eqre~^'^^)'{Hqre~^'^'^)' , so the time averaged 
energy flux is 

+ KK)/^ = i-KK + K:^;r)/(2c^)- (6.7) 

Thus the quantity (16. 6p has the physical interpretation as this time averaged energy flux, 
and it is then natural that its sum over all terminal edges should be the time averaged 
electromagnetic energy converted into heat in the circuit. 



7 A correspondence between electrical circuits and 
a subclass of electromagnetic circuits 

At fixed frequency, linear electrical circuits correspond to a subclass of EM-circuits, 
namely those where there are a sufficient number of magnetic elements and these all 
have ^qrs = 0. Let us consider, for simplicity, an n-terminal planar electrical network 
with terminal nodes at the vertices of a polygon and with the remainder of the circuit 
lying with the polygon. If Jqr is the complex current flowing from node q to node r and 
these nodes have complex voltages Vq and then we have 

Jqr = Yqr{Vq-V,), (7.1) 

where Yq^ is the complex admittance (having non-negative real part) of the circuit element 
joining these two nodes. 

For example, one may consider the four terminal network of figure [9]^a) which has 
two internal nodes. To build an associated EM-circuit, the first step is to triangulate the 
network by adding additional edges with zero admittance, as illustrated in[9]^b). To each 
triangle formed by this triangulation (not containing any nodes) with vertices g, r and s 
we assign a constant kq^-s = oo- In the limit as kq^s —>■ oo the equation (13.71) reduces to 

T; = r,'', = r;, Vqr + Vrs + Vsq = 0. (7.2) 



Following the ideas of Engheta, Salandrino, and Alu (2005) and Engheta (2007) we at- 



tach to each edge qr a dielectric cylinder with constant gqr = iYqr/oJ, which therefore 
will have non- negative imaginary part. [If the circuit element is a capacitor, then this 
will correspond to taking a value of the dielectric constant gq^. which is real and positive; 
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Figure 9: Construction of an EM-circuit corresponding to the planar electrical 4-terminal 
network (a). The first step in (b) is to triangulate the network, and place appropriately 
valued dielectric cylinders (not shown) along the edges, and magnetic triangles with 
A; = oo in each triangle. Then one adds a vertex below the network, and magnetic 
triangles with A; = cxo on the four triangular sides. The four new edges, marked by 
thicker lines, are the terminal edges of the EM-circuit. 
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if the circuit element is a resistor, then this will correspond to taking gqr with zero real 
part and positive imaginary part; if the circuit element is an inductor, then this will 
correspond to taking almost real and negative.] 
The equation (13.131) then becomes 

^ m m 

Jqr = YqrVqr, whcrC Jqr = — ^ T^r = ^qr^ (7-3) 

s=l s=l 

in which m = 1 or 2 is the number of triangles sharing the edge qr, and s indexes each 
of these triangles. 

We next introduce an additional node below the network, and for each pair t and 
u of neighboring terminal nodes around the polygon we construct the triangle tuO with 
constant ktuo = oo. As illustrated in figure [HI^c). This implies we have 

TL = Tt, = T,l V,^ + V^o + Vot = 0. (7.4) 

The edges uO, with u = 1,2, . . .n are taken as the terminal edges of the electrodynamic 
circuit, and no dielectric cylinders are attached to them. The second equations in (17. 2p 
and (17. 4p imply that we can assign a voltage Vg to each node such that 

Vgr = Vg-Vr, Vo = 0, Ko = K. (7.5) 

Thus (17.31) reduces to (17. ip . Also the first equations in (I7.2p and (I7.4p ensure that the 
total current is a sum of loop currents. Therefore Kirchoff 's law that the sum of currents 
flowing into a node equals the sum of currents flowing out of that node is automatically 
satisfied. Thus the standard electrical circuit equations are satisfied. 

Now the terminal edge variables Vuo, u = 1, 2, ... n, are the voltages at the terminal 
nodes of the electrical circuit. Also it is easy to see that the terminal edge variable 
Juo = AuQ/{iu) is the net current flowing out of the electrical circuit from node u to node 
0. Thus the map W/(ico') is the Dirichlet to Neumann map of the electrical circuit. 

If the electrical circuit is non-planar, then we modify the circuit so that all the terminal 
nodes are at the vertices of a (not necessarily convex) polygon lying on a plane below 
the circuit. Then the circuit above the plane is appropriately triangulated by adding 
additional nodes if necessary. A magnetic element with A; = oo is inserted in each 
triangle and appropriately valued dielectric cylinders are attached to the edges. Each 
pair of neighboring terminal nodes on the polygon are then attached with a magnetic 
element having = oo to an additional ground node situated below the plane. Each 
edge between a terminal node and the ground node is a terminal edge of the resulting 
EM-circuit. 
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8 Electromagnetic ladder networks, a characteriza- 
tion of their possible response matrices, and a ma- 
terial with non-Maxwellian macroscopic behavior 

Electromagnetic circuits can have many different topologies and seems very difficult to 
characterize their possible macroscopic matrices W, i.e. classify (for a given topology of 
terminal edge connections?) which matrices are realizable as the Schur complement of 
a matrix W° with elements (15.21) and (15. 3p . and which ones not. Here we restrict our 
attention to an important subclass of electromagnetic circuits, called electromagnetic 
ladder networks (EM ladder networks), for which such a characterization is possible. 
Consider, as illustrated in the simple EM-circuit consisting of two magnetic triangles qrs 
and rst joined by a cylinder with e = = along the internal edge rs. Assume they 
have the same constant kq^s = Kst = k. Then (13.71) implies 

= T;, = T,; = k{Vrs + Vst + V,r). (8.1) 

The edges sq and tr, labeled 1 and 2, are taken to be the terminal edges. They are 
without dielectric cylinders, so (I3.15P and (I3.13P imply 

^Iq = ^sq = —lu, = Atr = —hi, T!^s + ^rs = 0- (8-2) 

Suppose there are dielectric cylinders along the internal edges qr and st with the same 
constants Qqr = Qst = 9- At these edges (13.131) implies 

T^^ = uj'gVqr, T:, = u^gVst. (8.3) 

Solving these equations for I12 and I21 in terms of Vi = Vgq and V2 = Vtr gives 

h2 = -l21 = ki2{V2-V,), (8.4) 

where 

ki2 = l[l/k-l/iuj'9r' (8.5) 

has non-positive imaginary part, because k has non-positive imaginary part and g has 
non-negative imaginary part. From now on we ignore the internal edges of this simple 
EM-circuit, treating the simple EM-circuit itself as a basic ladder network element. 

The relation (18. 4p is similar to that associated with an element in an electrical circuit, 
although the physical interpretation of the variables is completely different. In the setting 



of an electrical circuit, using the notation of Milton and Seppecher (2008), 1 and 2 label 



two nodes, Vi and V2 are the potentials at these nodes, ilu/^ is the current flowing 
from node 1 to node 2, while ihi/oo is the current flowing from node 2 to node 1, and 
ki2 = 1/ L for an inductor, ki2 = —uJ^C for a capacitor, and k\2 = —iuj/R for a resistor, 
where L is the inductance, C the capacitance, and R the resistance. 
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Figure 10: A simple EM-circuit which is the basic element for constructing EM ladder 
networks. The triangles qrs and rst arc magnetic elements, and dielectric cylinders are 
attached to the edges qr and st. A cylinder, not shown, with e = /i = is attached to 
the edge rs to join the two magnetic elements. The edges rt and qs (which need not be 
coplanar) are terminal edges. 
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(a) (b) 

Figure 11: Figure (a) shows an EM ladder network which is the analog of the classical 
Wheatstone bridge, although its physical behavior is completely different. The two ter- 
minal edges are marked by the thicker lines, and for simplicity the dielectric cylinders 
along the edges of each basic ladder network element are not shown. Dielectric cylinders 
could be inserted along the edges between neighboring basic ladder network elements, 
in which case the EM ladder network becomes equivalent of a Wheatstone bridge with 
each node connected by a capacitor to ground. In figure (b) we have twisted the central 
bridge element so that it has a response governed by (18. 9p 

Building upon this analogy we can join a set of these simple EM circuits together, to 
obtain what we call an EM ladder network, as illustrated in figure [TlTa). An n-terminal 
EM ladder network consists of n + m edges Qa indexed by a = 1, 2, . . . , n + m, having 
no vertex in common. Each pair of edges {Qa-, Qfs) may have a simple EM-circuit (of the 
type just discussed) joining them, and from (18. 4p . we have the relation 

lap = = ka/siVp " K), (8.6) 

where kap is the constant associated with the simple EM-circuit, and ka/B = if there 
is no simple EM-circuit joining the edges. The first n edges Qa are the terminal edges 
of the EM ladder network (not to be confused with the terminal edges of the simple 
EM-circuits, which are all the edges Qa), and the remaining m edges are internal edges. 

Each pair of edges {Qa, Qp) may have an elementary EM-circuit of the type just 
discussed joining them. Each edge Qa, may have a dielectric cylinder, with constant Qa 
attached to it. If this edge is an internal edge of the EM ladder network then from (I3.13P 
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we have 

n+m 



E 



/a/3 = 9aVa, (8.7) 

/3=1 

(in which we set I^a = 0) while if this edge is a terminal edge of the EM ladder network 
then from fl3.15p we have, 

n+m 

Aa+Yl ^"/^ = -^^9aVa. (8.8) 

/3=1 

The response of the EM ladder network is then governed by the relation A = WV 
between the terminal variables A = (^41, A2, . . . , A^) which measure the free currents at 
these edges, and the variables V = (Vi, V2, . . . , Ki) which measure the tangential electric 
field at these edges. 

Equations fl8.6l) - fl8.8p are the same as those for electrical circuit in which the nodes 
may be connected to ground by a capacitor. It then follows directly from the results 



of Milton and Seppecher (2008) that for any fixed real frequency 00 any real symmetric 
matrix W may be realized by an EM ladder network having real positive values of the 
constants kajj and Qa, and any complex symmetric matrix W with ImW > can be 
realized by an EM ladder network having real positive values of the constants /c^/? and 
complex values of the constants Qa having non-negative real and imaginary parts. Thus 
at fixed frequency we have a complete characterization of the possible response matrices 
W of EM ladder networks, both in the lossless case, and in the lossy case. We also 
have a complete characterization of the possible response matrices W associated with 
the class of electromagnetic circuits where no two terminal edges are connected, since an 
EM ladder network can be constructed having these edges as its terminals, and having 
the desired response matrix W. 

We can introduce another basic ladder network element of an EM ladder network. 
Just by reversing the roles of the vertices t and r in the original basic ladder network 
element and using (13.111) . we obtain a basic ladder network element with the response 

Il2 = l21 = -kl2iVi + V2). (8.9) 

Of course utilizing such basic ladder network elements, as done in the example of figure 
[TlT b). does not enlarge the class of possible response matrices W of EM ladder networks 
at fixed frequency (which is already as large as possible without the introduction of such 
elements). 

Let us now sketch how one could get a material with non-Maxwellian macroscopic 
behavior using electromagnetic ladder networks. In the same way that one can build a 
cubic network of resistors so too can one build a cubic EM ladder network of basic network 
elements with the response fl8.6l) joined at edges Qa, with no dielectric cylinders attached 
to these edges (so that all ga = and it corresponds to the cubic network of resistors). 
Just as the cubic network of resistors responds macroscopically as a material with some 
effective conductivity, so too will the cubic EM ladder network respond macroscopically 
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in a way which is definitely non-Maxwellian. Without being too specific, for a periodic 
ladder network one will have some relation of the form K = SW, where V(x) is a 
suitably scaled local average of the variables Fq, S is the matrix governing the effective 
response, and i^^j(x), i = 1,2,3, is a suitably scaled local average of the variables Iaf3 
taken over the subset of basic ladder network elements which are "aligned" parallel to the 
Xj-axis. If such a cubic EM ladder network is embedded in a large cube having n = oo 
and e = with the terminal edges exposed at the boundary of the cube, then the interface 
conditions between the electromagnetic fields outside the cube, and the fields K(x) and 
V{x.) inside the cube will presumably depend on the geometric microconfiguration of the 
terminal edges of the EM ladder network at the cube faces. Obviously there is much to 
explore here. 
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